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O ' Abstract 

Qs , A system of nontopological solitons interacting through scalar and vector meson ex- 

change is used to model dense nuclear matter. The models studied are of the Friedberg- 
Lee type, which exhibit dynamical bag formation due to the coupling of quarks to a 
scalar composite gluon field a. It is shown that the high density behavior of such 
^ I models in the Wigner-Seitz approximation depends essentially on the leading power of 

fi ■ the quark-cr coupling vertex. By insisting that the parameters of any soliton model be 

chosen to reproduce single nucleon properties, this high-density behavior then selects 
>v>( I a promising class of models that better fit the empirical results — the chiral chromodi- 

V^ ' electric models. It is shown that one can adjust the model to obtain saturation as well 

as an increase of the proton charge radius with nuclear density. These two phenomena 
depend on the presence of a scalar meson and were not found in other nuclear matter 
calculations based on Friedberg-Lee type models. 
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1 Introduction 



Ever since the advent of quantum chroniodynamics (QCD) it has been popular to describe 
the nucleon in terms of bag or sohton models. There are many versions of such models, 
characterized by two extremes: the MIT bag model j^, where the nucleon consists of just 
constituent quarks arbitrarily restricted in a given volume, and the Skyrme model [3, where 
there are no quarks and the nucleon is instead a topological soliton of the pion field. Other 
models interpolate between these two extremes, seeking to combine the obvious emphases on 
the structure of nucleons within the MIT bag model and on nuclear interactions via meson 
exchange within the Skyrme model. For example, chiral bag models surround an MIT- 
like bag of quarks with a Skyrme-like cloud of pions (and perhaps other mesons as well). 
Clearly, each such model attempts to balance between low-energy degrees of freedom — the 
mesons — and high-energy degrees of freedom — the quarks (and possibly gluons as well). 
Our aim here is to develop a model that as simply as possible, but without losing essential 
physics, combines quark and meson degrees of freedom. This task is made more difficult by 
the "Cheshire cat principle" 0, an extrapolation of results found from chiral bag models, 
which find that low-energy nucleon properties are largely insensitive to the size of the quark 
bag. As the bag shrinks, the meson cloud forms more and more of the nuclear structure. To 
select a model, we must study nuclear properties that distinguish between a bag of quarks 
and a cloud of mesons. 

An obvious testing ground for any bag or soliton model is dense baryon matter, where 
the structure of individual nucleons will become as important as the interactions between 
neighboring nucleons. Preferably, the model will have a dynamical formation of the bag or 
soliton, for then one can treat the transition to a quark-gluon plasma consistently. Surpris- 
ingly enough, the Skyrme model, which may be thought of as a chiral bag model with the 
quark bag shrunk to zero size, does predict a dynamical transition to a phase of solitons of 
fractional baryon number at high densities ||^; however, identifying this phase with a quark- 
gluon plasma is certainly problematic. Thus in this paper we study dense nuclear matter 
in a set of models known as non-topological soliton (NTS) models. In particular, we study 
the Friedberg-Lee (FL) model and a class of related models, the chiral chromodielectric 
(xCD) models 0, §], as well as extensions of these models (or, one might argue, approxi- 
mations to these models) that explicitly include couplings to mesons. These models, based 
upon general arguments from QCD, are characterized by the coupling of quarks to a scalar 
field a that has a nonzero vacuum expectation field. This field is understood to be a com- 
posite gluon field. The interaction between the quarks and scalar field leads to a dynamical 
confinement mechanism, with the quarks carving a hole in the background scalar field. The 
structure of this bag depends precisely on how the quarks couple to the scalar field, and 
it is this coupling that distinguishes the various models we study here. All these models 
reproduce single nucleon properties reasonably well — indeed, chiral bag models show that 
single nucleon properties are relatively insensitive to the structure of the quark bag. We 
must look at high densities, where the bags begin to overlap, to see the differences between 
models with different quark-cr couplings. 



In the last few years there also has been considerable interest in the application of the 
quark-meson coupling (QMC) model to study the nuclear matter equation of state, as well 
as medium effects on the nucleon structure and nucleon-meson coupling, motivated by the 
apparent success of the Walecka QHD models 0, |10| in describing the properties of nuclear 



matter. The QMC model, initially suggested by Guichon |TT|, consists of non-overlapping 
nucleon (MIT) bags that interact through the exchange of scalar and vector mesons in 
the mean-field approximation (MFA). This model has been generalized to include Fermi 



motion and center of mass corrections [O] and was applied to nuclear matter (see [13] and 



references therein) and also to finite nuclei [|I^, |T3|. Of course, the assumption that the 
nucleons can be regarded as non-overlapping bags is only valid at low density, where these 
models seem to capture the essential physics. Already at nuclear saturation density, however, 
the internucleon separation is comparable to the nucleon radius, and at higher densities the 
assumption that the bags do not overlap clearly breaks down. This assumption becomes 
even more questionable when a modification of the bag constant in nuclear matter is taken 
into account. Jin and Jennings |]l6l and Miiller and Jennings [^ have recently shown that 



the introduction of a density-dependent bag constant can reproduce the EMC effect and 
reconcile the QMC results with those of the Walecka QHD-I model . As a result |T^ the 



bag radius grows with increasing density and the overlapping of the bags starts just above 
nuclear saturation density. 

To take into account the effects from overlapping nucleon bags — that is, to study a 
nuclear liquid rather a nuclear gas — it is clearly of interest to introduce some dynamics 
in the confining mechanism. Thus we are led to replace the MIT bag model by a NTS 
model. Although the study of nuclear matter properties was begun long ago (see P] and 
references therein) and was carried out with different soliton models and with different levels 



of sophistication (see |T^, IS] and references therein), none of these calculations included the 
effect of background meson fields on nuclear matter. Now, in principle, the non-topological 
soliton models — in particular, the xCD model, which in its full form includes perturbative 
gluon exchange — contain sea quarks and meson exchange explicitly. However, in practice it 
is very difficult to deal with anything besides the nuclear constituent quarks. In the interest 
of simplicity (and at the sacrifice of consistency), we extend our NTS models to include 
explicitly meson degrees of freedom. (Such an extension has been referred to as the local 
uniform approximation to xCD in ]|19| .) 

In the following, then, we shall study a system of non-topological solitons interacting via 
the exchange of scalar and vector mesons within the MFA. In modeling nuclear matter each 
soliton is centered in a spherical Wigner-Seitz (WS) cell. It has been argued that the choice 
of a spherical WS cell is more appropriate for a fiuid phase than a crystal as it represents an 
angular average over neighbouring sites. We consider two further approximations: in one, 
the WS calculation is simply used to provide an effective nucleon mass and the kinetic energy 
is then taken to be that of a Fermi gas, thus providing the correct low-density limit. A second 
approximation uses a Bloch-like boundary condition to calculate the band structure of the 
quark states. Neither of these approximations is quite satisfactory, and the second paper of 
this series is devoted to improving the modeling of a liquid of solitons. Nevertheless, we can 
still hope that the major qualitive features of dense nuclear matter will be reproduced even 
given the two approximations used here. Indeed, we find some encouraging results such as 



nuclear saturation and an increase of the proton rms with nuclear density. These effects are 
shown to derive from the background meson fields, which have a considerable effect on the 
formation of energy bands in nuclear matter. 

In this paper we limit our attention to Friedberg-Lee soliton type models extended to 
include the meson fields. The non-topological soliton models are presented in Sec. |[ The 
Wigner-Seitz approximation is then discussed in Sec. |. In Sec. § we study in some detail 
the trivial solution of the model. Using these solutions we show that the high density limit of 
the Friedberg-Lee type models depends on the leading power of the quark-a coupling vertex. 
The numerical results are presented in Sec. ^ 



2 Soliton models with quark-meson coupling 

Our starting point for studying the high density behavior of soliton matter is the Friedberg- 
Lee 0] non-topological soliton model, and we also consider related models like the chiral 
chromodielectric model of Fai, Perry and Wilets 0. In their simplest versions, these models 
include only constituent quarks and a single scalar field o that couples to the quarks. For 
the light quarks we shall assume tti^ = m^ = 0. Extending these models in the spirit of the 
quark-meson coupling model, we introduce in addition two meson fields: namely, a scalar 
meson and a vector meson V"^, which play important roles in quantum hadrodynamics. We 
assume these mesons couple linearly to the quarks. There is some freedom in the structure 
of the quark-meson vertex, in regard to its dependence on the soliton field a. Using the 
Nielsen-Patkos Lagrangian |2^, Banerjee and Tjon |^ have recently argued that in NTS 



models the quark-meson coupling should also depend on the scalar soliton field. Similar 
conclusions are reached by Krein, et al. [|l9l. However, we have found that this causes 
unwanted behavior within the mean field and Wigner-Seitz approximations used here (see 
Sec. 4), and so we report calculations that use a-independent quark- meson couplings. As 
we shall see, this choice has the advantage of reproducing the Quantum Hadrodynamics 
equation of state at low densities. 

Thus we take the Lagrangian density to have the form 

-r\d,ct>d^ct> - \ml<p' - \f,,F^^ + \mlV,V^ (1) 

where the a field self interaction is assumed to be 

f/(-) = 1-^ + ^-^ + ^^' + ^ ■ (2) 

The constants a, h and c are fixed so that U{(j) has a local minimum at a = (inflection 
point if a=0) and a global minimum at cr = cr^,, the vacuum value. The mass of the glueball 



excitation associated with the a field is given by racB = \/U"{a. 



The quark-o" coupling is taken to be 
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where the chroniodielectric function K,{a) has the form 

K{a) = 1 + 9{x)x''^[n^x - K + 1)] 



X 



a/ay 



(3) 
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In the following we will take n^ = 3. In perturbative calculations that include gluons the 
dielectric function K{a) is regularized in order to handle infinities in the one gluon exchange 
diagrams associated with a vanishing dielectric constant [0. Such a regulation is also useful 
for our numerical work. Koepf et. al. |2^ use the prescription 



K[a 



Kia){l 



+ Ky 



(5) 



where Ky is a constant often simply fixed at Ky = 0.1. However, this regulation forces 
g\ay) = 0, and we prefer instead to regulate as follows: 

K{a) = 1 + 9{x)x'^'''[ni^x — (n^ + 1 — Ky)]. (6) 

We find that properties of isolated solitons are independent of Ky for values even as large 
as .1, but results become sensitive to Ky at high densities. Here, we shall treat Ky as an 
additional parameter. 



The Euler-Lagrange equations corresponding to (|l|) are given by 

b^ {id^ + gvV^) - {nif + g{a) - g^^)] ^j = 

di,d^'a + U\a)+g'{a)ijij = 

d^dy + ml<P - g,^^ = 
-d.F^''-mlV'' + gy^Y^ = 



(7) 

(8) 

(9) 

(10) 

— — and g'{a) = -^^ . We solve Eqs. (f7|-|lO|) in the mean field approxi- 



where f/'(a) - ^^ c...^ y ,., - ^^ 

mation: we replace the soliton field a by a c-number a — 

their expectation values in the nuclear medium (p — >< (j) > 



a(r) and the meson fields by 



and Vii 



>< V^ >-- 
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with 00 and Vq constants. The approximation that the scalar and vector mesons fields can 
be regarded as constants while the soliton field is to depend on spatial coordinates stems 
from the long range nature of the light mesons and the short range nature of the soliton field 
due to the large glueball mass. In essence, the mesons are fast degrees of freedom and the 
glueball is slow, and we use a Born-Oppenheimer approximation. The resulting equations 
for the quark and the scalar soliton field are. 



-id-V + gyVo + j3 {mf + g{a) - 5f,0o)] ^fc = Cfc^fc 
-V^a + U'{a)+g\a) Y. V^fc^fc = • 

k(valance) 

We consider only valence quarks in our calculations. 



(11) 
(12) 



3 The Wigner-Seitz approximation 



In order to use Eqs. (pAl - P^ ) for the study of high density nuclear matter, it is sufficient to 
concentrate on a unit cell and solve these equations with the appropriate boundary condi- 
tions. In this paper we shall assume that each unit cell contains a single nucleon, that is, 
our unit cell coincides with a Wigner-Seitz cell. The single nucleon energy Ej:^ is a sum of 
two terms, the energy of the n^ = 3 quarks and the energy carried by the scalar soliton field 
a: 



En = riqeg + dr 

JWS cell 



\{Vaf + U{a) 



(13) 



The quark energy e^ should be regarded as the eigenvalue of Eq. (pJ] ) for isolated bags, or 
as the average energy in the band for dense nuclear matter. It should be noted that E^ 
cannot be identified with the nucleon mass as it contains spurious center of mass motion |p. 
In order to correct the center of mass motion in the Wigner-Seitz cell the nucleon mass at 
rest is taken to be 



Mn = ^El- < PI >ws , (14) 

where < P^^ >ws= ^q < p'^ >ws +^gb < (Va)^ >ws- The notation <>ws stands for an 
average over the Wigner-Seitz cell. Thus < p^ >ws is the expectation value of the quark 
momentum squared and nicB < (Va)^ >ws is the scalar soliton momentum squared. The 
latter average is obtained using a coherent state approximation [Q, which is essentially a 
single mode correction to the classical soliton mass. We find this latter correction to be 
small with respect to the total mass of the soliton and ignore it henceforth (it is not clear 
how valid this single-mode approximation is: at higher density we find that this correction 
exceeds the energy of the a field, with both going to zero). 

At low density the band width vanishes and the quarks are confined in separate bags. 
Then we expect the following approximation to be most accurate: we assume the individual 



nucleons move around as a gas of fermions with effective mass M^ given by Eq. (0), and 
so we get the following estimate for the total energy density at nuclear density p^: 

^ = ^r dkk'.jMj, + k^ + \ml<Pl - \mlVl (15) 

where 7 = 4 is the spin-isospin degeneracy of the nucleons. The Fermi momentum of the 
nucleons is related to the baryon density through the relation 

PB = ^4 • (16) 

The total energy per baryon is given by Eb = £ / Pb- The constant scalar meson field 0o is 
determined by the thermodynamic demand of minimizing £, which gives 

7 f^F o T&iEN- < PL >) 



4^'"^' ^° jEl~<P^^>ws+k'' 



The vector meson field is determined by averaging the Euler-Lagrange equation, Eq. (10) 
on a Wigner-Seitz cell yielding 
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These equations are similar to those of quantum hadro dynamics, the difference here being 
that the nucleon now has structure and thus the meson fields couple to the nucleon through 
its quarks. At low density the nucleon mass approaches its free value, and our mean field 
equations reduce to those of quantum hadrodynamics. 

The approximation used here for the nuclear system has often been adopted in modeling 
soliton matter, and is generally known as the Wigner-Seitz approximation. Each soliton 
is enclosed in a sphere of radius R such that ^R^ = 1/pb- On a periodic lattice the 

quark functions should satisfy Bloch's theorem iplf+a) = e^'''""ip{f) for the lattice vectors a. 
Concentrating on a single cell the Bloch theorem gives boundary conditions for the quark 
spinors in that cell. Although one can solve these boundary condition in a self consistent 
manner [|18|, we shall make the simplifying assumption of identifying the bottom of the lowest 
band by the demand that the derivative of the upper component of the Dirac function 
disappears at R, and the top of that band by the demand that the value of the upper 
component is zero at R |jT^. The quark spinor in the lowest band is assumed to be an 
s-state 

and the resulting Euler-Lagrange equations for the spinor components are 

duk 



(19) 



dr 



+ [rrif + g{a) - 5f,0o + {^k + 9vVo)] Vk = Q 



dvk 2vk . , . 

-r + + b"/ + 9{(^) - 9s 

dr r 



- (cfe + 9vVo)] Uk = . 
The corresponding equation, (|12|) , for the soliton field assumes the form 

U'ia)+g'ia)p,{r) = 0. 



-VV 



The quark density Pg and the quark scalar density ps are given by 
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(20) 
(21) 

(22) 

(23) 
(24) 



where the band is filled up to k. The quark functions are normalized so that there are 
three quarks in the Wigner-Seitz cell. The boundary conditions for the soliton field are 
cr'(O) = cr'(-R) = 0. The boundary conditions for the quark functions at the origin are given 
by u{0) = Uq and v{0) = 0, where Uq is determined by the normalization condition 



R 



Anr'^dr{u{rY + v{ry 



(25) 



The boundary conditions ai r = R are given by 

u',{R) = => VbiR) 
for the bottom of the lowest band, and 

ut{R) = 



(26) 



(27) 



for the top of the band. Using these equations we can solve for the corresponding e^ and e^. 
We assume the tight-binding dispersion relation 



e/c 



tb + [tt - eb) sin 




(2^ 



and that the band is filled right to the top kt. The assumption of such a dilute filling has been 
made previously ^j, and we do not discuss it further. The quark functions corresponding to 
the energy e^ can be simply obtained by integrating Eqs. ( pOD and ( pi] ) for each intermediate 
value efc. Substituting the dispersion relation into Eq. (p!3D, the nucleon energy is given by 



E 



N = -rj- I dkk'^ek + 



R 



Anr dr 



2 ^ ' 



+ U{a] 



(29) 



which is then used in (|T^) to determine the equation of state (p^S]). 



At higher density, where the bags begin to overlap, there is no reason to assume each 
quark is tightly bound to a single bag nor to impose that 3g groups move collectively with a 
well-defined momentum. Our equation of state cannot then be considered a good approxi- 
mation at higher densities, and we must then find a different approximation for handling the 
kinetic energy. This is the subject of the next paper in this series. For now we are content 
with studying low density behavior, with a special interest in whether our approximations 
can still produce saturation at the expected density. In fact, as one approaches the empiri- 
cal nuclear saturation density, the different models begin to distinguish themselves. This is 
discussed in detail in the next section. 



4 The Trivial solution 



The NTS models we are considering have a uniform plasma phase that is preferred at high 
densities. This corresponds to the solution cr = 0, so that the soliton bags "dissolve" and 
the quarks are free. For the original FL model, this solution is favored at unreasonably 
low densities. Moreover, in the WS calculations [§], for cell radii below ~0.8-0.9 fm only a 
trivial constant solution can be found. For the more sophisticated analysis reported in ||18|| , 
which includes higher partial waves in the quark wave functions and direct computation 
of the bands, the calculation breaks down already at a cell radius ^1.4 fm. We view this 
as a fault of the model and not of the WS approximation, for regardless of whether our 
restricting the soliton to a single spherical cell represents a dense system well or not, our 
model of the nucleon should nevertheless allow the soliton to be squeezed to volumes well 



below nuclear density before it breaks apart. (The saturation density of nuclear matter is 
Po=0.17fm~^, which corresponds to a WS cell radius i?o=l-12fm.) This leads us to look for 
improvements upon the original FL model, and the two particular extensions studied in this 
paper are generalizations of the quark-glueball coupling g{a) and the addition of explicit 
meson degrees of freedom. Before proceeding to the numerical results, let us look at how 
these extensions affect the trivial solution. 

The trivial solution to the WS equations (pO|-|2^) is Mfc(r) = uq, Vkir) = and cr(r) = o"o, 
where the constants Uq and ctq are independent of both r and k. Strictly, for k > this 
solution does not satisfy our boundary conditions; however, we find that below a given cell 
volume the numerical solution develops a singularity at the boundary in order to reproduce 
this preferred solution. Moreover, the more accurate self-consistent calculations of []TB| for 
the FL model (which do not assume just s-wave states in the quark wave function) show that 
the width of the band also narrows sharply at the onset of the trivial solution, indicating 
that the /c-dependence is not important. Thus in order to find the trivial solution that 
characterizes the breaking down of the WS calculation, we may assume that all the quarks 
are at the bottom of the band. This corresponds to squeezing an isolated soliton, and we 
insist that our model favors a nontrivial solution until the cell radius gets well below that 
corresponding to nuclear density. 

In the present section we shall consider altering the quark-meson couplings in our models 
as suggested in Lagrangian |^ ^, namely: 



9s 



gsg{(y), Qv -^ Qvoicr). 



(30) 



This will allow us to investigate whether the coupling to meson fields can help cure the 
breaking down of the FL model at too low density. 



(31) 



For the trivial solution, then, the normalization condition ( ^SQ gives 

3 xV2 



Uf) 



AttR^ 
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Pb ■ 



From Eq. (l2l|) we find the quark eigenenergy is 

e = mf + g{ao) [I - fi-.^o + 9vVo] , (32) 

where ctq obeys Eq. (j5^), modified to account for the change in q-a coupling, which leads to 

U'iao) = -3pBg'iao)[l - gs(po + 9vVo\. (33) 

Solving for the mean field values of 0o and Vq, we have 
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-^Pb9{(^o) and Vq = -^Pb9{(^o) 
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(34) 



yielding finally (for m/=0) 



e = 5'(o-o) { 1 - 3pB5f(o-o) 



mt 
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(35) 



and the total energy density 



^ = 3pB5'(oro) 



■^pB9{(ro) 



mt 



mt 



+ U{a,). 



(36) 



This is clearly a spurious solution corresponding to putting all the quarks in the lowest level. 
One needs to add the kinetic energy correctly to get the true energy of the quark plasma. 

We would like to select a model for which the trivial solution is found only at densities 
much higher than nuclear density, and we would like the quark mass e to be zero in the 
preferred phase. This latter condition implies o"o=0. However, for the FL model g{a) = g^^a, 
and therefore ao=0 is not a solution of Eq. (|33D . Using Eqs. (|3^ and substituting for U',g' 
and g, for the FL model Eq. (|33| ) becomes 



a^o + 2^0 + g^o 



' ' "' - -SpB9a 



o (9! 9l\ 

^PB9a \^^ ^0^0 

V mi mi / 



(37) 



If we turn off the meson fields gs,gv -^ 0, we see that the solution to (|37|) goes as —p^ at 



large densities. Not only does this blow up as p^ — ;> 00, but it gives the quarks an unphysical 
negative mass e ~ —gaPs ■ On the other hand, if we turn on the meson fields and make the 
usual choice of parameters so that gl/m^l > gl/m^l (which is necessary for saturation), there 
is now a positive high density solution that goes as p^^. However, there is still an unphysical 
negative solution, now diverging as pb- Thus we do not expect the inclusion of the meson 
fields to cure the problem of the FL model. 

Now let us consider a modified FL model in which we take the quark-glueball coupling 
to be g{a) = gaCr^- We shall call this the FL^ model. Then instead of Eq. (|3^) we have 



acTo +2^' 



O + g^O 
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For this model the trivial solution ctq = exists in the Wigner-Seitz approximation. There are 
also nonzero solutions, however. With the meson fields switched off there are two nontrivial 
solutions 



ctq 




8c , 

— {a + Gpsga) 



(39) 



which exist only at densities pb < a(36^/8ac — l)/6(y'o-. For the parameter choice used here, 
this corresponds to cell radii R > L67 fm. These solutions give the quarks a positive mass 
and do not necessarily signal any problems with the model. For nonzero q-a coupling, the 
solutions exist also at high density, behaving as ctq oc ±p^ . The quark effective mass (|35D 



vanishes in this limit. Thus the FL^ does not seem to have the problems of the FL model 
when using the WS approximation. 

Finally, let us consider the xCD models. From Eqs. (^ and (^ we see that for small 
a the coupling g{a) is of leading order ra^ in a. This means that, regardless of whether the 
mesons are present, for the xCD model with n^ > 2 there will be a solution to Eq. (p3D 
such that (To = 0. This corresponds to the desired free massless quark plasma phase, with 



10 



vanishing meson field averages, since then g{0) = g'{0) = 0. Furthermore, one can show 
that as pb ^ oo the only other solution has ctq ~ p^^ and gives the quarks a positive mass 
that also vanishes in the high density limit. Thus we expect the x^D models to provide a 
more reasonable description of dense nuclear matter than does the FL model, for there is no 
unphysical trivial solution to the WS equations that will cause problems. 

Indeed, the type of behavior predicted by our analysis of the trivial solution is reflected in 
our Wigner-Seitz calculations. As we show in Fig. 1, the xCD model exhibits Wigner-Seitz 
solutions down to very small cell radii (we never reached a breakdown point, taking R as 
low as O.OSfm), whereas the FL model breaks down at the expected point R ~ 0.8fm. The 
FL^ model can be taken down to lower cell radii than the FL, breaking down at i? ~ .25fm, 
well above the transition to the uniform plasma phase. However, the behavior of e^ below 
R = Ifm for the FL^ model does not seem desirable. This could be cured by allowing the 
quark-meson coupling to depend on a, but we find then non-smooth transitions to new non- 
trivial solutions in all the NTS models discussed here. We find more reasonable results if we 
keep the quark-meson coupling independent of a, and thus the x^D seems the preferable 
model. We shall present results only for this latter model in the next section. 



5 Results 

We have performed the calculations using a straightforward numerical integration of the 
equations of motion for the q and a fields. The mean field 0o is found by directly minimizing 
the free energy ([T5|). As a check, we have also performed some calculations with a relaxation 
routine, but this is far more time consuming. The numerical calculations with the xCD 
model were carried out using the following set of parameters 

a = 50fm-2 , b = -1300 fm-^ , c = 10^ g^ = 2, k^ = 0.1 . (40) 

The parameters of the potential U{a) are chosen to give a reasonable bag constant 
B = 46.6MeV/fm^ and glueball mass mcB = 1.82GeV. The value of a^ is calculated to 
be .285fm~^, the single soliton energy is 1391 MeV and the nucleon mass is 1176 MeV. 
The nucleon rms radius is 0.876 fm. The parameters were chosen to fit the rms radius and 
nuclear matter properties, resulting in a somewhat high value of the nucleon mass. We did 
not spend much effort in fine-tuning the parameters, however. As the short range repulsion 
provided by the vector field in QHD is provided by the Wigner-Seitz boundary conditions 
in the liquid soliton model, in what follows we shall take the coupling constant between the 
vector meson and the quarks to be zero, g^ = 0. On the other hand we shall treat gs as a 
free parameter in order to study its effect on the soliton matter. 

For completeness, we also list the parameter sets used for the other models studied in the 
previous section. For the FL model, we take the parameters as in Birse et. al. [0. (This 
set was also used in Ref. ||TB[|.) The full set of parameters for the soliton model is. 



0.0,6= -700.43 fm-i , c = 10\^<, = 10.98 . (41) 
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These parameters correspond to a single soliton energy 1260 MeV, and nucleon mass 902 
MeV after removing the center of mass motion. For the FL2 model we use the same potential 
parameters and a coupling g^j = 60fm~^. 

We proceed to our presentation of the results for the xCD model. To our knowledge, 
this is the first application of this particular model to the study of dense matter. Previous 
work has dealt with few nucleon systems f^. 

In Figs. 2 and 3 we show the fields u{r) and a{r) for several values of the cell radius R. 
As can be seen from its behavior at the cell boundary, the quark field Uk{r) displayed in Fig. 
2 is that corresponding to the bottom of the band k = 0. Fig. 3 shows that the depth of the 
bag decreases only slighty as the cell radius is lowered, indicating that the quarks remain 
tightly bound. (Note that for the FL model, however, the bag depth actually increases with 
smaller R [p!7| §.) That the quarks are tightly bound in the bag can be seen clearly in Fig. 
4, where we display the quark density Pq{r) given by Eq. (|23|). There is only a small overlap 
with quarks from neighboring cells even at radii below 1 fm. In Figs. 1-4, the quark-meson 
coupling has been set to zero. 

In Fig. 5 the top and bottom of the lowest energy band is shown as a function of R 
for several values of the quark-meson coupling Qg. The point at which the band begins to 
form, R ^ 1.6fm, is not very sensitive to the value of Qs- On the other hand, the structure 
of the band depends strongly on the coupling to the scalar meson, becoming wider as Qs 
is increased. In Fig. 6 we present the energy of the soliton E^ given by Eq. (^9]), and in 
Fig. 7 we show the total energy per baryon Eb of the system, derived from Eq. (^). As 
expected, the introduction of the scalar meson to the the soliton matter results in attraction 
and saturation. The empirical nuclear saturation density corresponds to Rq = 1.12 fm, 
whereas from Fig. 7 we see, for example, that our NTS nuclear matter energy per baryon 
has a minimum at i? ~ 1.35fm for gg = 1. The saturation energy is Eb — M)^^ ~ 20 MeV, as 
compared to the empirical value Egat = 16 MeV. The compression modulus of nuclear matter 
is K = R^d^Es/dR^ at the equilibrium point: for the present model we find K ^ 1170MeV. 
This is to be compared to empirical estimates that usually lie in the range 100-500MeV, with 



K = 200MeV the generally accepted value [ 23 1. These results should be regarded as only 



order of magnitude, as at higher density we are surely underestimating the kinetic energy 
of the system. We shall develop a more reasonable model of the liquid state, which leads to 
a more quantitatively accurate EOS for solitonic nuclear matter, in the next paper of this 
series. 

The solution for the xCD model is stable down to very low values of R (at least as low 
as 0.05 fm), and it is apparent that with proper calibration this model can be used as an 
excellent starting point for the study of nuclear matter. Another important feature of this 
model is the increase of the nucleon rms with density. Fig. 8. This unexpected outcome of 
the model is in accord with the EMC effect. Checking the sensitivity of this effect to model 
parameters, we found that the increase in the rms is insensitive to recoil corrections, but 
might disappear with different choices of the soliton parameters. 



12 



6 Discussion 

In this paper we have investigated a class of non-topological sohton models that generalize the 
Friedberg-Lee model. We have studied nuclear matter in the Wigner-Seitz approximation, 
where neighboring bags begin to overlap at higher density, as a means of distinguishing 
between the various models. The models differ in the precise form of the coupling between 
the quarks and the scalar gluon field, and in the presence and form of couplings to explicit 
meson fields. Of the models studied, we have found the chiral chromodielectric model to 
exhibit behavior best in line with phenomenological expectations. For this model, there is no 
breakdown in the Wigner-Seitz calculation as the cell radius is decreased, thus overcoming 
a serious shortcoming of the original FL model. We consider it a computational necessity to 
consider only constituent quarks and ignore colored gluons and sea quarks, and thus these 
restrictions along with the addition of explicit meson fields may be seen as an approximation 
to the original theory. The original x^D model is closely based upon QCD, exhibiting 
absolute confinement, and it is therefore satisfying to see that this model is selected by our 
studies of nuclear matter. 

We include mesons along the lines of quantum hadrodynamics, and we use the mean 
field approximation. In our final calculations we have considered only the scalar meson, but 
it is straightforward to include also the vector meson, which should be done for fine-tuning 
the parameters. We note that we found it best to use a quark-meson coupling that is not 
modulated by the quark-glueball coupling, in contrast to that argued for in Refs. |]IU|, |2T]]. 
Indeed, when using a g — coupling proportional to g{cr), we find jumps to qualitatively 
different solutions as density is increased. The resulting curve of energy as a function of cell 
radius has several bumps, which is clearly an undesirable feature. However, this should not 
necessarily be taken as an argument against models that employ quark-meson couplings that 
depend on a, but rather merely a statement that when using the mean field approximation it 
is more consistent to use a a- independent q-(j) coupling Qg. Indeed, one can view this coupling 
as a "mean field approximation" to a more fundamental coupling — namely, Qs ~ {gsg{o'))- 

In the end, then, we have used a xCD model with cr-independent coupling between 
quarks and the scalar meson. We find that the inclusion of the scalar meson provides a clear 
saturation point for nuclear matter and that a rough fit to both empirical nuclear matter 
and single nucleon properties can be obtained. One of the most interesting features found is 
an increase in the nucleon rms radius at intermediate densities, in line with the EMC effect. 
This is dependent on the presence of the scalar meson. There are clearly several refinements 
to our model and our calculations that must be considered before fine-tuning the parameters. 
Among these are the inclusion of other mesons in the model, the addition of perturbative 
gluonic effects, and an improved calculation of the quark wave function along the lines of Ref. 



IS]. One can also improve our handling of the spurious CM motion. Foremost, however, we 



need to improve our modeling of the liquid state. 

In this paper we have used a low-density approximation in modeling nuclear matter 
with solitons. This consists in using a Wigner-Seitz approximation to calculate an effective 
nucleon mass. Then the kinetic energy is added to the system by taking the motion of the 
nucleons to be that of a Fermi gas. The assumptions of this picture include: the nuclear 
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medium restricts any given nuclear bag to a spherical cell; any given nucleon moves slowly, 
so that it can be constructed at rest and then boosted; the quarks remain tightly bound 
inside the bag; the nucleons move independently within the medium. Only the first of these 
assumptions can remain valid at high densities, where the nuclear medium will form a liquid 
in which each nucleon's motion is localized on short time scales and quarks need not remain 
tightly bound inside the bag. Since we are ultimately interested in studying the transition 
to a uniform quark plasma within our non-topological soliton model, we must do a better 
job of modeling the liquid state at high densities. This is the subject of our next paper. 
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Figure 1: Energy of the bottom level of the band as a function of cell radius for the 
Friedberg-Lee (FL), modified Friedberg-Lee (FL^) and Chiral Chromodielectric (xCDM) 
models 
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Figure 2: Upper component u{r) of quark wave function at the bottom of the band for 
several cell radii R in xCD model 
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Figure 3: Soliton bag a{r) for several cell radii R in xCD model 
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Figure 4: Quark density Pq{r) for several values of R in xCD model 
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Figure 5: Energy of the bottom and top levels of the band versus R for several values of 
the quark-meson coupling Qs 



1500.0 



1000.0 



> 



W 
T3 

03 



500.0 



0.0 




E, gs=o 
E, gs=o 
Eb gs=i 
E, gs=i 
E, gs=2 
E, gs=2 



0.0 



1.0 



2.0 



R[fm] 



3.0 



20 



Figure 6: Soliton energy En (before subtracting spurious CM motion) for several values of 
9s 
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Figure 7: Energy per baryon Eb of nuclear matter as function of R and Qs in the xCD 
model 
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Figure 8: Nucleoli rms radius as function of R and Qs in the xCD model 
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